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We used the Hartree-Fock approximation to classify the electronic phases that might occur in 
a transition metal nanowire. The important features of this situation are orbital degeneracy (or 
near-degeneracy) and interactions favoring locally high-spin configurations. In this circumstance, 
spin density wave and triplet superconductivity states are favored. If the interactions favor locally 
low-spin configurations as in the previously studied spin ladder systems, orbital density wave and 
singlet superconductivity are observed. 
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I. INTRODUCTION 

Self-assembly epitaxial techniques have enabled the 
fabrication of one-dimensional (ID) atomic wires com- 
posed of adatoms confined at step edges on surfaces of 
appropriately chosen substratea^"^ In many physically 
relevant cases, the surface bandgap structure of the sub- 
strate material is such that the electronic states of the 
adatoms are decoupled from the bulk substrate bands 
(at least to leading order) and a one-dimensional electron 
system can be realized. STM (Scanning Tunneling Mi- 
croscope) and ARPES (Angle-Resolved Photoemission) 
measurements of Au nanowires grown on the Si(577) sur- 
face show that a charge-density wave (CDW) occurs at 
low temperatures^!. The physics of Au nanowires is still 
the subject of debate^', but it appears that the rele- 
vant band is derived from Au s-states for which electron- 
electron interactions are relatively weak, and the domi- 
nant physics may be associated with lattice instabilities 7 . 
However, many other adatom/substrate pairs are possi- 
ble, and this opens new possibilities including the study 
of one-dimensional electron gases formed from transition 
metal d-orbitals. For example, wires compo sed of Co 
adatoms have been furrowed on a Cu substrate^. 

From the theoretical point of view, the important fea- 
tures of transition metal-based wires are the orbital de- 
generacy of the transition-metal d-levels, which permits a 
rich set of on-site interactions and the small size of the or- 
bitals, which leads to larger interaction effects. In partic- 
ular, the Hund coupling favors locally high spin configu- 
rations, potentially leading to interesting spin structures. 
Motivated by these ideas and the recent experimental 
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success^, in this paper, we investigate the physics of a 
nanowire in which the important electronic states are de- 
rived from the transition metal d-orbitals. We use mean- 
field theory to establish the phase diagram and elucidate 
the general classes of behavior. While mean-field theory 
is not an exact description of interacting electron systems 
especially in low spatial dimensions, it should tell reason- 
able indication of what physics is relevant and provide a 
starting point for more exact treatment. A subsequent 
paper will use renormalization group and bosonization 
methods to obtain a detailed picture of the same model 



when one-dimensional nature is significant. 

Consideration of the physics of transition metals leads 
to models with multiple electronic bands with more or 
less arbitrary interactions. Models of this general class 
have been previously considered in the literature, both 
for their intrinsic interest!^, and as steps toward under- 
standing heavy fermiorP^lEI] and high temperature su- 
perconductor systems^H^. In these models the multiple 
bands arise from physically different atoms: in the heavy 
fermion case, one band represents the local moments and 
the other the wide band of conduction electrons; in the 
high T c case an important motivation has been models 
of "spin-ladder" compounds^. While a formal mapping 
may be established between these models and the mod- 
els of interest here, the different structure of the interac- 
tions leads to different physical behaviors. Models more 
directly analogous to those of present interest have also 
been investigated, and the focus has been put on ferro- 
magnetism and orbital ordering in the strong coupling 
regime^HIII xhe relation of these works to the results 
obtained here is discussed below in Sections [TXJ |IV| and 

YD 

The organization of this paper is as follows. In Sec. [i^ 
we explain the model and the symmetries. In Sec. |III[ we 
define order parameters encoding the physics arising from 
the backscattering and forward scattering. The general 
properties of mean-field solution and approximation we 
employed will be presented in Sec. |IV| In section |Vj we 
show the obtained phase diagrams in the weak-coupling 
limit 



and discuss the results. Sec 
diagrams obtained in the strong coupling limit 
Sec. |VII| is conclusion and summary. 



VI is devoted to phase 
Finally, 



II. MODEL 

While many of our considerations are general, we 
are motivated in particular by recent success gro wing 
monatomic Co chains at step edges on a Cu surface^. 
Band calculations indicate that the Co electrons are hy- 
bridized into the bulk Cu states mostly away from Fermi 
energy, and Co electrons form one-dimensional partly 
filled bands near Fermi energy^. Thus, as a first approx- 
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(b) Atomic wire with d-orbitals 



FIG. 1. Schematic pictures of two-leg Hubbard ladders and 
d-orbital model. In two-leg Hubbard model, there is hopping 
along the ladder, t, and transverse hopping, t±. The on-site 
interaction is U. For atomic wire with d-orbitals, we have in- 
traband hopping, t.4 and tB, with various on-site interactions, 
U, U', and J. 



imation, we may consider that the system is described by 
a multi-orbital Hubbard-like model representing the Co 
d-orbitals with local onsite Coulomb interactions 

H = E E - l T' (4 ms c jm > s + H.c.) + H int . (1) 

m,s 

Here c^ s is the annihilation (creation) operator for a d- 
electron in orbital m with spin s at site i. t™ m is the 
hopping between from orbital m on site i to orbital mf on 
site j. The interaction terms Hi nt will be discussed below. 
Through the paper, we set the lattice constant equal to 1. 
The presence of the surface breaks the symmetry between 
d levels and may lead to an arbitrary ionization level. 
For the sake of simplicity, we will consider here only the 
case where Fermi energy crosses two orbit als, m = A, B, 
although in the general case one would have 5 d-derived 
bands with an arbitrary fermi energy. Furthermore, the 
rotational symmetry in Hi nt as we will see always allows 
us to diagonalize the hopping matrix, so we will ignore 
t AB . 

In the weak-coupling limit, the band structure is char- 
acterized by four Fermi points: two Fermi momentum, 
kA and kg, and two chiralities, r = R,L. R{L) rep- 
resents electrons around positive (negative) Fermi mo- 
menta. The total particle number is n = 2(kA+kB)/^- It 
is useful to distinguish the two cases of half filling (n = 2) 
and arbitrary filling (n ^ 2). For arbitrary filling we ex- 
pect in general that ^ but two special cases 
are possible: we may have degenerate bands, Ha = &b; 
alternatively, one of the two bands may be half- filled i.e. 
kA = § 7^ kg. If n — 2 then fc^ + kg — ir; in general we 
expect kA ^ fcs but both bands could be commensurate, 



kA — kg — § ■ When two bands have equal Fermi mo- 
mentum and Fermi velocities, the kinetic term acquires 
0(3) ~ SU{2) orbital symmetry. 

For the two-orbital system the interaction terms have 
the following form: 

-Hint = U E n imt n iml 
i,m 

+ V E n iAsniBs 
i,s 

+ (V - J) E n iA s riiBs (2) 
i,s 

— J E c iAs C iAsc\ Bs CiBs 

+ Jl E ( c Lt c tAi. c uH ciB t + H - c -) 

i 

where and rii ms — c\ ms Ci ms is the electron density and 
s = — s. U and U' indicates on-site Coulomb repulsion 
between two electron in the same band or different bands, 
and J represents Hund coupling favoring high spin state. 
J' is the so-called pair-hopping term. For a transition 
metal ion in free space all of these parameters are pos- 
itive. We assume that the symmetry breaking by sub- 
strate primarily affects the hopping terms in the Hamilto- 
nian without changing the local orbitals too much. This 
enables us to use free-space rotation symmetries to re- 
duce the number of interaction constantP^. In this case 
we have 

J = J' (3) 
U = U' + 2 J. (4) 

The first equality is derived from the fact that Wannier 
wave functions are real, and the second one represents 
rotational invariance in orbital space. With this sim- 
plification, the interaction terms now have U(l) orbital 
rotational symmetry about y-axis. For a transition metal 
ion in free space, U ^S> J > 0, so that all interaction pa- 
rameters are positive. Screening will reduce the value of 
U, but will lead to only negligible changes in j^22l, and 
most calculations indicate that even the reduced value of 
U is greater than J. 

Now we talk about the symmetry of the Hamiltonian 
in Eq. and Eq. @. The model has SU(2) symmetry 
about rotation of spins. The U(l) rotational symmetry 
in orbital space in interaction terms is preserved only 
when two bands are equivalent. Otherwise, the orbital 
symmetry of the total Hamiltonian is completely broken. 

Lastly we compare our models to previously stud- 
ied ones. In particular, this model is similar to two- 
leg Hubbard ladd er models, which have been studied 
extensiveljEmSEIIEll, The comparison is sketched in Fig. 
[I] In essence the two sides of a rung of the ladder (or, 
more precisely, the odd and even parity combinations 
of these states when there is strong transverse hopping) 
map to the two atomic states we consider. Generically, 
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in the ladder problem, the hopping across a rung is non- 
vanishing, implying in our language fc^ ^ ks- However, 
the interactions of two models are quite different. In 
two-leg ladder problems, the intra-chain interactions U 
and transverse hopping t± are supposed to induce effec- 
tive antiferromagnetic coupling (—U/t\ < 0) b etween 
two sites connected by a rung of the laddePEStH. On 
the other hand, the model considered here doesn't have 
such hopping since two states in an atom are orthogonal, 
but instead it has Hund coupling (J > 0), which favors 
high-spin states. 

More recently several groups attempted a general clas- 
sification of the physics of ladder systems with generic 
interactions either with transverse hoppingf- 19 ! 20 !, or with- 
out transverse hopping^. When two Fermi momenta are 
different, ^ k B , our model is considered in Refs. IT§1 
and Effil while when fc^ = fcg, we have pair- hopping term 
which is not included in Ref. 1371 These studies employ 
perturbative renormalization group and bosonization as- 
suming weak interactions and one-dimensionality. While 
we will pursue this direction in subsequent paper, here 
in this paper we will use a mean-field approach, which 
allows us to access the strong coupling regime of the 
model - which is hard to success by perturbative RG and 
bosonization - as well as to classify the different possible 
states. 

More directly related studies has been done on the 
same mode l focusing on ferromagnetism and orbital 
ordei 



In the strong coupling regime, U ^> J >■ t, 
both analytical and numerical calculation show ferro- 
magnetism and orbital antiferromagnetism as a ground 
state at quarter-filling. Around quarter-filling, ferromag- 
netic state is robust to hole-doping and electron-doping 
with nearest neighbor hopping. However, below quarter- 
filling, the inclusion of further hopping leads to disap- 
pearance of ferromagnetism while the state still exists 
above quarter-filling. Haldane gapped state with 5 = 1 
is expected at half-filling. 



III. ORDER PARAMETERS 

Mean-field theory involves minimizing the energy with 
respect to a free-fermion density matrix characterized by 
non-vanishing expectation values for fermion bilinears in 
the particle-hole (ph) and particle-particle (pp) channels. 
A general bilinears is characterized by a momentum, q, 
spin, s, and orbital, to. We distinguish between q = 
cases (ferromagnetism, ferro-orbital order, and supercon- 
ductivity), and q ^ (charge/spin density wave, orbital 
density wave, FFLO or pair-density wave superconduct- 
ing states). We will mostly not be interested in the de- 
pendence of the expectation values on the magnitude of 
the fermion momentum and will mostly be interested in 
electronic states near the fermi points. Therefore, we la- 
bel the bilinears by the chirality, spin, and orbital indices 
without explicitly denoting the q or fermion momentum. 



The basic objects are: 

/A \ss' \mm' 1 \ ^ / + \ 

\<-±ph) rr i — / J \ c rms c r'm' s' / 



and particle-particle bilinears, 



/ a \ss':mm 1 \ * /ft 

(A PP ) rr , = jj 2^ (mscl ms cl — 



(5) 



(6) 



where Cm S is the annihilation (creation) operator of elec- 
tron with chirality r, orbital to, and spin s. We will use 
the following convenient basis to represent these, 
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Y, ^^(A.^-'—'+H.c. (7) 

mm' ss' 

^pp 

mm' ss' 



r mm ,a 3 ss ,(A ph r ;mm (8) 



where i,j = (0,1,2,3) and r and a are Pauli matri- 
ces with r" b = a {) ab = S a b. These operators trans- 
forms as tensor of rank 2 under the rotation of SO (4) ~ 
SU (2) spin x SU (2) OI . bitaJ . Due to the SU (2) spin symmetry 
of the Hamiltonian, we can take quantization axis along 
z-direction for spins, and consider only j = and 3. 
However, we will keep all 2 = ~ 3, since there are cases 
without any orbital symmetry. er°( 3 ) combination gives 
spin singlet (triplet), that is, charge (spin) mode. Sim- 
ilarly t°( 3 ) gives orbital singlet (triplet), and t 1 ^ gives 
symmetric (anti-symmetric) combinations of orbitals. 

We first discuss the cases with r = r' in Eq. ([5| , which 
corresponds to spatially uniform density order. We con- 
sider only ferromagnetism (FM), orbital- ferromagnetism 
(OFM), and combinations of these two. These orders are 
characterized by non-zero density polarizations: 



&fm — 



"At) 



n B t) - \n 



Aqfm = (rut) + ( u m) ~ ( n Bf) - (nBi) ^ °- 



(9) 



In Sec. |VIJ we will consider FM(+OFM) state, and this 
is a state where Afm is maximized first, and then A.ofm 
is maximized next. 

Next, we turn to the density wave for r ^ r' . In or- 
der to classify these phases, we label our particle-hole 
order parameters by the phases at each Fermi point and 
transferred momentum; there arc 4 possible cases for in- 
terband (q — kA + fcs) and intraband (q = 2k m ) or- 
der (Table |H| . s-wave has the same phases at all Fermi 
points, p^-wave changes its sign under parity transfor- 
mation, R -f-> L, and p^-wave does under band exchange, 
A <-> B. d-wave is odd under both transformations. Ap- 
plying this classification, we find that i = and I are 
both s-wave, although the former is intraband type and 
the latter is interband type. We put " ' " for interband 
order to distinguish these two. i = 2 and 3 is found to 
be interband and intraband p y -wave accordingly. 

When a band is commensurate, we have another family 
of order parameters called "bond" order (BOW), which is 
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(i,3) 


Particle-hole order 


Particle-particle order 


P 


L 


S 


(0,0) 


Charge density wave (CDW) 


d'-wave singlet SC(d'SS) 


-1 








(0,3) 


Spin density wave (SDW) 


p' y -wa,ve triplet SC (p' y TS) 


1 





1 


(1,0) 


s'-wave charge density wave (s'CDW) 


Pji-wave singlet SC (p y SS) 


1 


1 





(1,3) 


s'-wave spin density wave (s'SDW) 


d-wave triplet SC (dTS) 


-1 


1 


1 


(2,0) 


p'y-wave charge density wave (p'CDW) 


s-wave singlet SC (sSS) 


1 


1 





(2,3) 


Py-wave spin density wave (p'SDW) 


Pa;-wave triplet SC (p^TS) 


-1 


1 


1 


(3,0) 


p^-wave charge density wave (pCDW) 


s'-wave singlet SC (s'SS) 


1 


1 





(3,3) 


Py-wave spin density wave (pSDW) 


p^,-wave triplet SC (p^TS) 


-1 


1 


1 


TABLE I. Classification of order parameters. " ' " 


indicates that the order is interband type. 


The 


eigenvalues 


of each 



superconducting phase under parity (P), orbital rotation (L), and spin rotation (S), are also listed. Here the orbital (i) and 
spin (j) indices are defined in Eq. ^ and Eq. Q. Particle-hole channels are even under parity. 



basically the density- wave slid from on-site to "on-bond", 
and is the same as dimcrization. The only difference 
between site order and bond order is phase of the order 
parameter; A is real for on site order, and imaginary for 
bond order. We found that the energy gain is maximized 
when order parameter is real, indicating that always on- 
site order has lower energy. Therefore, we will ignore the 
bond orderings in the remainder of the paper. 

To be complete, in two-leg ladder problems, p'CDW 
is more commonly called as orbital antiferromagnet 
(OAF^m, or staggered flux (SF) stated, pydensity 
waves can be called orbital-density waves, but should be 
distinguished from PDW in Ref. [19j which is a bond- 
order. 

We label the particle-particle channels for supercon- 
ducting orders in the same manner (Table [I]). In partic- 
ular, when kA 7^ fcs, the order parameter with i = 0, 3 
has non-zero momentum, q = kA — ks] the order exhibits 
periodic structure in real space s imilar to FFLO (Fulde- 
Ferrell-Larkin-Ovchinnikov) stat d 40 * 4 ^. However, in our 
case, there is no external field to split the spin up and 
down electrons. This possibility of FFLO state in multi- 
orbital system without external field was first pointed 
out by Padilha and ContinentincS^. The d'SS state often 
appears in two-leg ladder problems^. 



IV. METHOD 

We employ the standard Hartree-Fock approxima- 
tion, reducing the quartic part of the Hamiltonian to 



"Angular momentum" 


(A,R) 


(A,L) 


(B,R) 


(B,L) 


s 


+ 


+ 


+ 


+ 


Px 


+ 




+ 




Py 


+ 


+ 






d 


+ 






+ 



TABLE II. Angular momentum and phase at each Fermi 
point. 



quadratic, AB ^ (A)B + (B)A - (A)(B), where (A) and 
(B) are determined by minimizing the energy. These 
expectation values correspond to (quasi) long-ranged 
orders^ induced either by forward scattering or by 
backscattering. Since it is not feasible to treat all the 
scattering processes in whole parameter space, we mainly 
investigated two regimes: the weak-coupling regime 
where back-scattering is dominant, and strong-coupling 
regime where both forward scattering and backscattering 
compete. This separation is motivated from the obser- 
vation that Stoner's scenario of phase transition driven 
by forward scattering requires coupling to be larger than 
some critical value, although the backscattering always 
opens a gap even in the weak-coupling limit in ID. There- 
fore, in the weak-coupling regime, we ignore forward 
scattering, and focus on back scattering. In the strong 
coupling regime, we first assume that forward scattering 
drives the system to some kind of density polarization, 
and then consider the effect of residual backscattering to 
these polarized states. 

A. Effect of backscattering 

In this subsection, we explain the treatment of 
backscattering since we can treat it in the same man- 
ner for both weak and strong coupling regime by us- 
ing constant density of states (DoS). Although the de- 
tailed form of DoS is important to determinate the phase 
boundary between strong-coupling phases and weak- 
coupling phases, this approximation is justified within 
each regime: in the strong coupling regime, kinetic terms 
are less important than interactions; in the weak-coupling 
regime, electrons far away from Fermi energy is irrele- 
vant. 

We first focus on a single band case. The quadratic 
Hamiltonian obtained by mean-field approximation can 
be diagonalized, and the system is gapped at Fermi en- 
ergy, cp. Using the energy e measured from tp, the new 
dispersion is found to be, 

±v/e 2 + .g 2 A 2 (10) 
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where A is order parameter in Eq. ([5| and Eq. 
and g is the corresponding coupling constant. In our 
model, g is expressed by some linear combinations of U 
and J. Complete list of coupling constants expressed by 
U and J is given in Table III of the Appendix. Note that 
there exist contributions from Umklapp processes at half- 
filling, and extra interband scattering when k^ = kg. 
Under the assumption of constant DoS, the energy gain 
by this gap is given by 



SE 



J (e+ Ve 2 + ff 2 A 2 ) de - gA 2 



(11) 



where v is the density of states for single band, and A is 
the cut-off or bottom of the band. The values of these 
parameters are different in the weak-coupling and strong- 
coupling regime, and we will explain it below. 

By maximizing the energy gain in terms of A, we get 
the analytical solution to the gap equation, 



A 
.9 



sinh 



(12) 



The stable condition for ordered phase is g > 0. In order 
to obtain the phase diagram, we compare the energy of 
possible phases, and choose the order which gives small- 
est energy as the ground state. Thus, phase boundaries 
indicates 1st order transition from one minima to another 
without coexistent region. 

When multiple bands are involved, the calculation be- 
comes more tedious. In our two band model, the im- 
portant possibilities are: two bands have the same or- 
der independently by intraband scattering, or two bands 
have an order by interband scattering. For the former 
case, the Eq. (11) remains the same where v is density 
of states for each band, and total energy only depends 
the averaged density of states, v intra = (i/^ + z/g)/2. For 
the latter case, the dispersion becomes more complicated 
in general, though the final result depends only on sin- 



gle parameter, 



' 'inter 



v B l )/2. Thus, there is 



an inequality between interband and intraband density 
of states, 



^intra ^ ^inte 



(13) 



Therefore different density of states lead to suppression 
of interband scattering. 

Next, we see the weak-coupling limit and the strong 
coupling limit of the above results. 



and the energy gain for single band by gap-opening is 
found to be 



SE ~ vg 2 A 2 



(15) 



The density of states is fixed to be the value at Fermi 
energy, and the cut-off A is taken to be small compared 
to band width, At. The density of states at Fermi energy 
is connected to the Fermi velocity of the band, v, and 
it is given by v = (2ttv) 1 . Therefore, we see that the 
velocity difference suppresses interband processes from 
Eq. " ' 



(13). 



C. Strong coupling limit 

In the strong coupling regime, along with the backscat- 
tering, we have forward scattering which induces static 
orders or density polarized states. We assume that this 
polarization is maximized in the strong coupling limit, 
and study the effect of back scattering on each polarized 
state. We will not consider partially polarized states, 
which might appear between non-polarized state and 
fully polarized state, because the possible intermediate 
phases are complicated and depend sensibly on the de- 
tails. 

As a constant density of states, we will use the averaged 
value for tight binding dispersion, v = l/(4t), since most 
of the electrons participate in density-wave formation in 
the strong coupling limit. The energy of each state con- 
sists of two parts: static energy, and energy reduced by 
backscattering. The former is simply given by the sum of 
kinetic terms, and static density-density interaction. As 
we take constant DoS to be l/(4t), kinetic term becomes 
2tn(n — 1) with particle density n in each band. The re- 
duction of energy by density wave formation is obtained 
from Eq. (JTTJ) . In particular, in the limit of g — > oo, it 
becomes 



SE 



1 2 a 2 1 A 2 

-qu A vK 

4 y 2 



with 



A~ vK 



3v 2 g 2 



(16) 



(17) 



The backscattering tries to use all the electrons to form 
a density wave in the strong coupling limit, so the cut-off 
A in Eq. ( 11 ) is taken as the energy of band bottom. 



B. Weak coupling limit 



to 



In the weak-coupling limit [g 0), Eq. ( 12 ) is reduced 
2A 2 

A^ — e"^, (14) 
9 



V. WEAK-COUPLING PHASE DIAGRAMS 

Here we present Hartree-Fock phase diagrams in the 
weak-coupling regime. In order to obtain the phase di- 
agrams, we compare the energy of possible phases, and 
choose the one with lowest energy as ground state. Along 
with the order parameter with (quasi) long-range corre- 
lation, these phases are characterized by the number of 
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FIG. 2. Hartree-Fock phase diagrams for away from half-filling. Physically relevant region is U > and U/2 > J > 0. The 
number in the parenthesis indicates the number of massless modes in charge (C) and spin (S) sectors. We used v = 1/4 to get 
the phase diagram for Ua 7^ fcg = tv/2. 



gapless excitations in charge and spin modes. We denote 
a system with m massless charge modes and n massless 
spin modes as CmSri^. Without any interaction, the 
original Hamiltonian has 4 bands and this corresponds 
to C2S2. 

We first explain the three cases away from half-filling, 
and then see the phase diagrams for systems at half- 
filling. For all the cases, the physically relevant parame- 
ter region is U 3> J > 0, although we investigated various 
parameter regions beyond this restriction. 



A. Away from half-filing 

The weak-coupling phase diagrams where filling is 
away from half-filling are given in Fig. [2] for three cases. 
We will look at each parameter region, and explain the 
dominant physics which governs the phase. 



1. U > J > 

First we concentrate on upper right plane (U, J > 0) 
since this is the physically relevant parameter region. In 
the small J region, the Coulomb repulsion U dominates 
the physics, and as in the one orbital Hubbard modeP^ 
the ground state is a spin density wave. In the generic 
case of two incommensurate fermi wave vectors the phase 
of the spin density wave is not pinned between the two 
channels and there is a continuous family of solutions. 
In kj± — fcs case, the two Fermi momenta are the same, 
and the relative phase mode is pinned down; there is no 
degeneracy here. 

In terms of number of massless modes, total charge 
mode and relative charge mode are both massless in de- 
generate SDW phase where both spin modes are mas- 
sive. Thus, the degenerate SDW phase is expressed as 
C2SQ. This represents two independent metallic spin- 
gapped chains of (7150. Non-degenerate SDW phase in 



k_A = ks case has massive relative charge mode so it 
becomes C1S0. 



2. J>|(/|>0 

The condition for Coulomb integrals to be positive 
(U 3> J > 0) means that this regime is unlikely to be 
relevant to real materials. One notable feature for large 
J is that we have p-wave supercon ducti vity, which is also 
observed in numerical calculations^^. This is different 
from the case of two-leg ladder systems, where purely 
repulsive Coulomb interaction leads to d-wave supercon- 
ducting ground state of spin sing letP, the p-wave super- 
conductivity is triggered by attractive interaction. We 
can understand this by looking at the limit of J — > +00, 
where the spin on the same site is fully polarized, but 
an orbital degeneracy remains. So the only on-site inter- 
action with dynamical consequences is (U — iJ)nAcrn-Ba 
in Eq. By employing the knowledge that negative- 
U Hubbard model has spin singlet superconductivity as 
the ground state^, we find that the analogous ground 
state of this limit is interband orbital singlet supercon- 
ductivity with parallel spin, which is namely p^TS. The 
degeneracy of p^TS with p^TS arises from the absence of 
pinning effect between two SC as is the case for SDW and 
pSDW. When two Fermi momentum arc not equal, these 
superconductivities show periodic modulation of order 
parameters in real space similar to that found in FFLO 
state. 

At last, the degenerate p^TS and p^TS state is C2SQ, 
and non-degenerate p' TS for %a = &s is C150. 



3. -(/>|J|>0 

Next, we turn our attention to negative U region. This 
regime is also not likely to be realized in transition met- 
als. When \U\ is large enough, it is naturally expected 
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(1, 0) 

-i.o< u — i.o r u 

-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

(a)k A ^ k B {h)k A = k B = tt/2 

FIG. 3. Hartree-Fock phase diagrams at half-filling. Physically relevant region is U > and U/2 > J > 0. The number in 
the parenthesis indicates the number of massless modes in charge (C) and spin sector (S). The sSS states should be read as 
S-Mott or S'-Mott state. 



that attractive U < gives some kind of superconduc- 
tivity; indeed, we found p y SS for positive J, and sSS for 
negative J when two Fermi momenta are different. p^SS 
is replaced to CDW when kA — ks- 

To understand these phases, here we consider two- 
particle local eigenstates. There are 4C2 = 6 locally 
possible states. The spin triplet (S — 1) states are 



\S = 1, S z = 1) = c J+ c 



10} 



\S 
\S 



o> 



1 



V2 



a ss' c As C Bs 



<|0) 



(18) 



-1) 



cti 1 c 



|0) 



The on-site energy is E$ = i = U — 3 J . Among three spin 
triplet (S = 0) states, £7(l) or bitai doublet states are 



1 



1 



|5 = 0, 2 ) 



(19) 



^2 cr mm' c mt t 'm'^ 



a|0) 



with Es=q,- = U— J. The last piece is <7(l) or bitai singlet 



15 = 0,- 



J J 



alo) 



(20) 



with Es=o,+ = U + J. This indicates that, for large 
negative U, interband superconductivity with S = 1 is 
preferable for J > 0, and intraband spin singlet super- 
conductivity is preferable for J < 0. The latter super- 
conductivity is indeed sSS phase in negative J region. 
On the other hand, the positive J region doesn't match 
withj)j,SS states in the phase diagram. This discrepancy 
is attributed to the different numbers of allowed scatter- 
ing processes; when two Fermi momenta are different, 
the number of interband scattering process is fewer than 
that of intraband ones. Therefore, interband ordering is 
suppressed. For example, the following interband process 
is prohibited when kA 7^ ks, 



nAariBs 



c \l<7 C ARo Cb Rs CBLa 



H.c. 



(21) 



although similar intraband process is allowed, 



"ma ' v m<r 



^mL<j C ™-R<y^ m fi s CmL(j 



H.c. 



(22) 



Therefore, \S = 0, z) is more suitable in positive J, and 
this corresponds to p y SS. Of course, when J becomes 
sufficiently strong, the energy gain by spin alignment be- 
comes predominant, and the system exhibits spin-triplet 
superconductivity. Similarly, the CDW phase in upper 
left area for equivalent bands shows up since it is strongly 
enhanced due to the "nesting" of fc^ = ks although Pj,SS 
is not affected. 

Here all the phases are non-degenerate, so only the 
total charge mode is massless, CI SO. 



4. -J>|[/|>0 

Finally, the large negative J region is again described 
by the nAansa term in Eq. ([2]), and ground state should 
be interband orbital density-wave, which corresponds to 
SDW of simple Hubbard model. So the possible candi- 
dates are either s'CDW, or pCDW. Taking into account 
the ordering of the fermionic operators in nAansa, we 
find that s'CDW has correct sign to be the ground state. 
For kA = ks case, s'CDW and pCDW are degenerate as 
^(l)orbitai symmetry requires. 



B. At half-filling 

The phase diagrams for half-filling cases are shown 
in Fig. [3| At half-filling, the most of the argument of 
the general filling still apply, but we have to take Umk- 
lapp processes into consideration. Since Umklapp pro- 
cesses enhances only density-wave states, superconduct- 
ing states which appear in negative U region are now re- 
placed by CDW as in kA = ks case away from half- filling. 
An interesting new phase is s'SDW which is located be- 
tween s'CDW and SDW for fc^ ^ case. At this spe- 
cial filling, interband Umklapp process is enhanced, so 



Py SS 






' SDW 




. pSDW 


^SS 









1.0 








0.5 




CDW 




0.0 






SDW 


-0.5 




sSS 




U -1.0 

- 


.0 


-0.5 0.0 





(a)fej^ ^ fcg ^ tt/2 at n ^ 2 (b)fcyl 7^ at n = 2 

FIG. 4. Phase diagrams for va/vb ~ 40. Physically relevant 
region is 17/2 > J > 0. 



s'SDW is dominant at small J. However, s'SDW is sta- 
ble only when U > J, although SDW is stable for all 
U, J > region (See Table [TTT| of the Appendix). Thus 



at large J, SDW is again dominant, and we obtain the 
above phase diagram. 

The sSS phase at half-filling should be read as S-Mott 
or S'-Mott state, which often appear in two-leg ladder 
problems^; at commensurate filling, we know that the 
system is insulating due to Umklapp process and not 
metallic. These Mott insulating states have similar order 
parameter as sSS except total charge mode when it is 
written in bosonic fields, and turn into sSS upon doping. 

Finally for most of the phases appearing at half-filling 
is completely gapped, (705*0, except a region where 
s'CDW and pCDW are degenerate. In this degenerate 
region, the orbital sector is massless, C1S0. 



C. Effect of velocity difference 



As we pointed out in Sec. |IV| the velocity difference 
suppresses interband scattering process, and intraband 
order becomes dominant. In our cases, the dominant 
phases appearing arep y SS, CDW, SDW, pSDW and sSS. 
As the velocity difference gets larger from va/vb ~ 1, 
the phases governed by interband scattering are gradu- 
ally excluded, and beyond va(b)/vb(A) ~ 40, which cor- 
responds to Vmtra/vinter ~ 10, the whole phase diagram 
is covered by intraband type orderin g (Fi g. |4|). The 
kA 7^ kg = tt/2 case looks like Fig. |4(a) when either 
va/vb or vb/va becomes large. Whe n kA = the 
phase diagrams are similar to Fig. |4(b)| regardless of the 
filling. 



VI. STRONG-COUPLING PHASE DIAGRAMS 
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FIG. 5. Schematic phase diagram for the strong coupling 
limit, 17 > J > f with J/t = 5 including inhomogeneous 
phases. Green: FM(+OFM), Orange: FM+ODW, Gray: 
mixed state of FM(+OFM) and FM+OAFM, Purple: mixed 
state of FM+ODW and FM+OAFM, Blue: mixed state of 
FM+ODW and AFM. 



(1) : FM, FM(+OFM), FM+ODW, FM+OAFM 

(2) : SDW, AFM, OFM+SDW, OFM+AFM, 

where AFM stands for antiferromagnetism, and OAFM 
is orbital-antiferromagnetism. In particular, we will dis- 
tinguish two spin-density waves: AFM with q = tt and 
SDW with q — 2k p. The two states are identical when 
kp = tt/2 but while the SDW state is driven by a fermi 
surface instability and is the only important state in the 
limit of weak coupling, the AFM state is stabilized by 
commcnsurability (Umklapp) effects and may exist for a 
range of carrier concentrations near the commensurate 
value. We similarly define orbital density wave (ODW) 
with wave vector q = 2kp, and OAFM with wave vector 
q = 7r. 

Calculations very similar to those leading to Eq. [TT] 
and Eq. [12] give the following results for q = ir orderings 
in a single band with n < 1/2. The ground state energy 
is given by 



E G s = v 



-2t(l-2n) 



(-Ve 2 + .g' 2 A 2 ) de + g'A 2 
+ (static interaction energy), (23) 



We turn now to the phase diagrams in the strong cou- 
pling regime, C/> J>i. Since we are mostly interested 
in high-spin state, we compared the energies of following 
4 ferromagnetic states and 4 antiferromagnetic (or SDW) 
states: 



with v = l/(4i). We used g' — 2g to emphasize that the 
coupling constant is doubled at q = tt due to Umklapp 
process, while static energy from q = part is just with 
g. The kinetic energy for n > 1/2 can be obtained by 
particle-hole symmetry. The solution for gap equation is 
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l) 2 + 1 + 2(2n - l)cosh 

\ v 9' 
(24) 

This goes to ~ in the strong coupling limit. The solution 
exists only when the density is close to half-filling, n c < 
n < 1 — n c , with 



1 



> I 2 \ , f 2 

' cosh — — 1 — cosh 



(25) 



This becomes 1/2 at g —¥ 0, and goes to as g —> oo. 
Thus, in the intermediate coupling, the density wave with 
q = n is stable only around half-filling. 

For simplicity we assumed two degenerate bands with 
constant density of states. Particle-hole symmetry allows 
us to investigate only < n < 2. Comparing the ener- 
gies of the 8 states discussed above, we obtained ground 
state phase diagram, which is given in Fig. [5] Below 
quarter-filling, FM(+OFM) - where only single band is 
occupied - is dominant with large U; this configuration 
does not cost any interaction energy. The transition be- 
tween FM+ODW and FM(+OFM) below n = 1 can be 
understood by Stoner's scenario where orbital sector be- 
comes polarized above critical value, U c . The precise 
behavior of the phase boundary as n — > depends on 
details; for example, DoS of an isolated chain diverges 
at very small n leading to the smaller value of critical 
interaction strength. 

As the filling becomes closer to quarter-filling, 
FM+OAFM is found to be stable because it can use 
Umklapp processes to cancel static interaction energy, al- 
though this solution is unstable if too much holes or elec- 
trons are doped (See Eq. ( 25 ) ) . As we plot the energies of 



these states, we found that phase separated state exists 
below n = 1, which mixes FM+OAFM and FM(+OFM) 
for large U and FM+OAFM and FM+ODW for small 
U . At exactly quarter-filling, the system is homogeneous 
FM+OAFM state. 

Above quarter-filling, the FM+ODW is more sta- 
ble than FM(+OFM) and FM+OAFM; the energy of 
FM+ODW in the strong coupling regime is roughly 



E 



1 



(U -3J)n 2 



(26) 



FM+ODW below half-filling, and at half-filling, the sys- 
tem is totally occupied by AFM. 

Now, we'd like to compare our mean-field phase di- 
agram to the previously obtained results. At quarter 
filling (n — 1), Kugel and Khomskii^l, and Cyrot and 
Lyon-CaerPl found FM+OAFM as the ground state by 
strong-coupling e xpans ion, and this is confirmed by nu- 
merical calculatio n 1 26 ! 27 !. This result can be understood as 
follows: when spins are totally aligned, Fermi momentum 
are doubled, and effectively the system is at half-filling. 
Then we may regard orbital index as pseudo spins, and 
the system exhibits pseudo-spin density wave. 

Away from quarter-filling, the Umklapp process is 
killed so we expect OAFM is less dominant; indeed, 
Sakamoto et alW^ found FM+ODW with tight-binding 
DoS. They also found that adding far-neighbor hopping 
to get constant DoS replaces FM+ODW to paramag- 
netism (S — 0) in n < 1 though the system remains to be 
FM in n > 1. This is because FM is induced by double- 
exchange mechanism for electron-doped case, but for 
hole-doped case, it is driven by purely one-dimensional 
"spin-charge separation"^, which is fragile to perturba- 
tion of far-neighbor hopping. These observations do not 
contradict with our result above quarter-filling, though 
we have FM(+OFM) instead below quarter-filling. We 
think the FM(+OFM) state is actually more or less sim- 
ilar to the paramagnetic state without double occupancy 
in Ref. [STJ since both configurations do not cost any 
interaction energy below quarter-filling. The mean-field 
treatment picks up FM(+OFM) among other configura- 
tions which do not have doubly occupied sites. On the 
other hand, the ferromagnetism in Ref. [57| is induced 
by spin-charge separation, which is not a phenomenon 
captured by mean-field theory. Therefore, we conclude 
that the ferromagnetism of FM(+OFM) in FigM5] and 
the state seen in the numerical results of Ref. 12 7l have 
different origins. 

At half-filling, the system is claimed to be Haldane 
type where fully-polarized spin 1 on each site are antifer- 
romagnetically coupled by exchange interaction. Slightly 
below half-filling, a phase separation between Haldane 
phase and FM+ODW was founds, which agrees with 
our results for small U. 



VII. SUMMARY 



while the energy of FM(+OFM) and FM+OAFM is lin- 
ear in n, E = (U - 3J)(n - 1). Thus, FM+ODW is 
energetically preferable above n ~ 1.18. Again the tran- 
sition from FM+OAFM to FM+ODW is smeared by a 
phase separated state of these two. 

Near half-filling, n ~ 2, AFM state appears in weak 
U regime by the same reason for FM+OAFM to appear 
around quarter-filling. However, the kinetic part of AFM 
does not cancel the static part completely, and the resid- 
ual interaction makes this state unstable as U gets larger. 
AFM state forms an inhomogeneous mixed state with 



In this paper, we used mean-field theory to determine 
the phase diagrams of the two-band Hubbard model 
for a wide range of interactions including the intraband 
Coulomb repulsion U, interband Coulomb repulsion 
U'(= U — 2J), Hund coupling J, and pair-hopping J. 
For transition metal ions, we expect U 3> J > 0. 

First, we looked at the weak-coupling regime where 
back-scattering is dominant. For equal Fermi velocities, 
we observed five general features irrespective of band 
structure, each corresponds to the following parameter 
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regions: 

(1) U ^S> J > 0: This parameter regime is relevant to 

real materials. In this region, various SDW or- 
ders (SDW, s'SDW, and pSDW) are most dominant 
similarly to simple Hubbard model. When kA ^ ks 
with incommensurate filling, SDW and pSDW are 
degenerate since there is no phase pinning effect be- 
tween independent SDWs in two orbitals due to the 
incommensurability. When &a + ks = 7r, the inter- 
band Umklapp process enhances s'SDW in small J 
region, and becomes dominant. 

(2) J ^> \U\ > 0: The ground state is spin triplet super- 

conductivity (p'TS) for incommensurate filling, and 
CDW at half-filling. The former is driven by the at- 
tractive interaction by large J and the Cooper pair 
is formed by interband electrons. When ^ fc^, 
their order parameters have non-zero momentum 
such as FFLO state^l but without external field. 
The CDW at half-filing is induced by the strongly 
enhanced Umklapp process. 

(3) — U >• — J > 0: This is again described by single- 

band Hubbard physics, and sSS order develops. 

(4) — U 3> J > 0: When filling is commensurate or two 

Fermi momentum are equivalent, these conditions 
allow additional scattering processes for CDW, and 
this becomes the ground state. On the other hand, 
for kA 7^ ks and incommensurate filling, intraband 
pSS is dominant due to positive J and suppression 
of interband process. 

(5) -J > \U\ > 0: The ground state is interband CDW 

(a' CDW), which is an orbital analogue of SDW 
in single-band Hubbard model. When the system 
has [/(l)orbitai symmetry, pCDW is degenerate to 
s'CDW. 

The velocity difference reduces interband process, and 
intraband ordering becomes dominant relatively. As we 
increase the velocity difference from va/vb = 1, we ob- 
served that interband type ordering gradually expelled 
from the phase diagram, and above va/vb ~ 40, the 
phase diagram is completely covered by intraband type 
ordering. 

Second, we investigated the strong coupling regime of 
the model. We found ferromagnetism is almost always 



achieved, and various orbital orders are realized depend- 
ing on density and interaction. Around quarter-filling, 
FM+OAFM with q — tt modulation is stable. The region 
above quarter-filling is dominated by FM+ODW with 
q = 2kF , and there exists phase separation between these 
two phases. Below quarter-filing, FM+ODW is domi- 
nant for small U, and this is replaced by FM(+OFM) 
for larger U. The transition between these two phases 
by changing the magnitude of interaction can be un- 
derstood by Stoner's scenario. As density gets larger, 
these states pass through phase separated regime and be- 
come FM+OAFM at quarter-filling . Close to the half- 
filling, there is an inhomogeneous mixed state of AFM 
and FM+ODW for small U. 

For the Co/Cu system, it is located in U J > re- 
gion if we neglect the hybridization between the Co wires 
and Cu surface. Although the real system has more than 
two orbitals, we may expect the following results about 
the Co/Cu systems from our calculations. First, if the 
system is in the weak-coupling regime, it exhibits SDW, 
s'SDW, or degenerate SDW-pSDW state depending on 
the band structure. In bosonization scheme, these phases 
are replaced by spin-gapped phases or quasi-long range 
order phase. In strong-coupling regime, we may have 
FM(+OFM), FM+ODW, or FM+OAFM. If the filling 
is very close to half-filling and the U is not large, AFM 
is also possible. Of course the effect of hybridization or 
larger number of orbitals will introduce more complex 
physics into the system, but we will not pursue it here. 

In a subsequent paper, we will present a result by 
renormalization and bosonization taking quantum fluctu- 
ation into consideration. The low energy effective Hamil- 
tonian obtained after integrating out the high frequency 
modes could be different from the microscopic Hamilto- 
nian we considered here, and different ground state is 
expected to appear. Especially several Mott insulating 
phases which are dual phases of density-wave should be 
investigated. These order parameters are expressed by 
non-local "string" operators, and are not considered here. 
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Appendix A: Coupling constants 



In this Appendix, we list up the coupling constants of the phases we considered in this paper (Table III ) . The stable 
condition for the gap equation requires the coupling constant to be positive. For commensurate filling, Umklapp 
process doubles the number of possible scattering for density waves, though the superconducting order is not affected. 
When only one band is commensurate, the intraband scattering is enhanced for density-wave formation. 
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TABLE III. Coupling constant for each phase 
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